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In the present paper we consider the problem of a horizontal magnetic dipole situated 
in or out of a lossy magnetoplasma halfspace when the magnetostatic field is parallel to the 
axis of the dipole. The rigorous formulation of the problem is carried to the point where 
the determination of the pertinent boundary coefficients remains to be a straightforward 
but not a simple algebraic process. Due to the prohibitive algebraic complexity involved 
in the explicit finding of the boundary coefficients, a high frequency approximation is intro- 
duced and the approximate form of the boundary coefficients found. The field integrals 
are evaluated in the air region for the condition when the observation point is at a large 
distance from the source. The results are applied to finding the radiation pattern in the 
air due to a magnetic dipole situated near the lower edge of the ionosphere. It is found that 
the most important correction due to the earth's magnetic field occurs near the interface 
and parallel to the axis of the dipole. 

1. Introduction 

The problem of a horizontal electric and magnetic dipole in air and in the presence of a 
conducting, homogeneous, and isotropic halfspace was originally formulated and solved by 
Sommerfeld [1926]. Extending the work of Sommerfeld to a conducting halfspace, Moore 
[1951] and also Banos and Wesley [1953] solved the problem of the horizontal electric dipole 
wdiile Wait and Campbell [1953] and Wait [1959] solved the problem of a horizontal magnetic 
dipole. 

The problem of radiation in free space from sources in magnetoplasma halfspace was first 
considered by Barsukov [1959]. He obtained expressions for the Poynting vector in air for 
the cases when the magnetostatic field was normal to the boundary and an electric dipole was 
either normal or parallel to the boundary. He also presented a numerical example for a vertical 
dipole at the boundary and for the particular set of conditions considered he found that the 
plasma's anisotropy manifested itself most strongly in the directivity of the radiation pattern. 
A similar problem was also considered by Arbel [I960]. 

In the present paper we shall treat the case of a horizontal magnetic dipole situated either 
in or out of the magnetoplasma when the magnetostatic field is parallel to the dipole axis. The 
choice of a horizontal magnetic dipole as the source of the electromagnetic waves was motivated 
by the fact that such an antenna in the presence of an interface effectively produces both 
vertically and horizontally polarized waves; thus, both wave polarizations could be studied 
in a single problem. Moreover, since the electric field lines form circular loops about the axis 
of a magnetic dipole and the magnetostatic field is parallel to the same, then the alternating 
electric field is normal to the steady magnetic field and, thus, maximum interaction of the 
electromagnetic wave with the magnetostatic field could be expected. 

2. Dipole in Magnetoplasma 

2.1. Rigorous Formulation 

The geometry of the problem is showm in figure 1. The horizontal plane 2=0 coincides 
with the interface between the anisotropic homogeneous plasma and air. For conven- 
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Figure 1. Geometry of the problem of a magnetic 
dipole in magnetoplasma. 



ience we shall call the plasma medium (1) and the air medium (0) . The plasma, in addition to its 
anisotropy, may be lossy while the conductivity of the air is zero. Moreover, we assume that 
both media have the same magnetic inductive capacity of free space, pt . The factor causing 
the anisotropy of the plasma, the magnetostatic field H DCj is oriented in the positive ^-direction 
as well as the source of the electromagnetic waves, the magnetic dipole. 

The definition of the present boundary value problem implies the solution to Maxwell's 
equations subject to the usual boundary conditions at the interface and proper behavior at 
infinity. In the air region it is convenient to introduce the Hertzian vector potential of the 
magnetic type whereas in the plasma region it is found to be more convenient to work with the 
actual field components. 

For the purpose of this problem it will be assumed that the source of the electromagnetic 
waves consists of a small wire loop carrying an electric current Ie~~ i<at . When the loop is small 
enough, its electromagnetic effects can be adequately represented by its equivalent dipole 
moment 



m=I-S 



(1) 



where S is the surface vector of the area enclosed by the loop. We localize the source by 
writing for the magnetic current density 



J m =md(x)5(y)5(z+h). 



(2) 



In the magnetoplasma region the magnetic field will satisfy the inhomogeneous vector 
wave equation 



(-VXe- 1 VX+k 2 o )H l = -io,e md(x)d(y)d(z+h)l x 
where Tis the permittivity tensor given by [Tyras and Held, 1959] 
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where N is the number of free electrons per cubic meter, |e| = l. 6X10" 19 coulomb is the 
electronic charge, and ra e =9.1X10~ 31 kg is the electronic mass. 

In the air region we use the Hertzian vector of the magnetic type defined by 



£-ia) Mo vxn (M) 

and satisfying the homogeneous vector wave equation 



(6) 



(v 2 +^)n W) =o. 



(7) 



The formulation of the present boundary value problem can be simplified a great deal 
by expressing the field components in the magnetoplasma and the air regions in terms of their 
triple Fourier integral representation in Cartesian coordinates in the transform space as well 
as in the configuration space. To this end we introduce a triple Fourier transform pair 
defined by 

F(a l ,a 2 ,a,)=~fr 2 fjJF(x,y, z)e' i( ^ + ^ + ^dxdydz 



F(x, y, 2) = ^yr 2 J j j F(a u a 2) a 3 )e i( ^ x+ ^ v+ ^ z) da 1 da 2 <la B . 



(8) 



In what follows we shall also need the transforms of the derivatives. These can be obtained 
by integrating by parts where the vanishing of the integrated part is assured providing that 
the fields have a proper behavior at infinity. 

Tims transforming (3) one obtains for the plasma region 
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(9) 



(10) 



The system of algebraic equations in (9) can be solved using Cramer's rule. The 
determinant of the coefficients of the square matrix on the left is found to be 



A=^(al-sV(alsl) 



where 
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(11) 



-2al±^/[(x-f)^- X ^ ^ ojJ+4J"«»*wV (12) 



The result in (12) is not surprising. It is typical of what could be expected from analysis 
of a double refracting medium like magnetoplasma. We have two modes of propagation 
characterized by the eigenvalues Si and ,s* 2 with which terms "ordinary" and "extraordinary" 
are often associated. These eigenvalues are given by the zeroes of A and are the poles of the 

3 The correspondence between our notation for plasma parameters and that of the magneto-ionic theory [Rat cliff, 1959] is as follows: e=l+C£7, 
v = CY t , f=l — X(U. We feel that our notation is more convenient for our purpose. 
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field components. We note, for instance, that if a 2 = 0:3=0 in (11), which is equivalent to 
saying that the wave does not vary along the y and 2 coordinates, we obtain for the eigenvalues 



(ai)i.2=*oV€±^ (13) 

which can be recognized as the propagation constants of the right- and left-hand circularly 
polarized plane waves propagating in the direction of the magnetostatic field. Alternately, 
if ai = a^=0 we find 

(a 2 ) 1 =k -y/x 

(a 2 ) 2 =£oVf (14) 

which again can be recognized as plane wave propagation constants transverse to the steady 
magnetic field. Setting ai = a 2 =0 and solving for o 3 gives the same results as in (14). 

We now proceed with formally solving the system of algebraic equations in (9). Using 
Cramer's rule and the results of (11) and (12) one obtains 
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where 

$i(a u a 2 ,a 3 ) = (xk 2 o— a?)[f(x&o— a 2 ) — x(a|+a|)] — k 2 $ol\ 

$2(01,02, o 3 ) = fx^o— f (1 — ^ 2 )oI— x(o|+Os). (16) 

The inversion with respect to the o 3 -transform variable can be performed immediately by 
integrating in the complex a 3 -plane along a path that runs along the real axis then along a 
circular arc in the upper half-plane for 2+A>0, and the lower half-plane for 2+A<0. It can 
be shown that the contribution from integration along the circular arc vanishes providing 
lm{s li2 }>0 and the resulting integral is equal to the sum of the residues at the poles s x and s 2 . 
Furthermore, inverting with respect to the o r and o 2 -transform variable one obtains the 
integral representation of the primary excitation as follows: 
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To satisfy the boundary conditions of the problem, we shall also need an appropriate 
complementary solution of the homogeneous system of (9). From the theory of linear differ- 
ential equations [Ince, 1956] it follows that the magnetic field satisfies 

-> 
(al-sl)(.al-sl)H M = (18) 
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where a z in this equation corresponds to the differential operator id z . The solution to (18) 
can be written immediately 
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where we discarded solutions with positive exponentials for physical reasons. The coefficients 
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dj are not all independent. For since H (c) must satisfy the given system the coefficients must 
be related by 
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where the square matrix elements F ti {si t ^ correspond to the square matrix of the original 
system (9) with a z =—s u or a 3 =— s 2 . Now we express all coefficients C ti in terms C n and 
C12 and normalize obtaining 
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Inverting with respect to the ai and a 2 transform variables we obtain the complementary 

solution as follows: 
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The coefficients A x and A 2 , representing waves reflected from the interface, will be determined 
uniquely from the boundary conditions. 

The Hertz vector potential in the air region must satisfy the wave equation (7). The 
appropriate solution can be written 
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n < M)= _W_ CC Bie Ha lX + a2 v+s z> daida2 (24a ) 

co 

n^> = "~" e C ((B a e"i*+''* +, <Pda 1 da i (24b) 

CO 

where s =(kl— ol\— al) V2 and we shall require that lm{s }>0. The coefficients B x and B 2r 
representing waves transmitted through the interface, will be determined uniquely from the 
boundary conditions. 

The boundary conditions to be satisfied by the Cartesian components of the field vectors 
require continuity of the tangential components of the electric and magnetic fields at the 
interface z=0. This implies the following 

k 2 U^ + d x (d x U x M) + D z Tl ( z M) ) =H xl 
tk 2 d y Ui M) = ^H 2l -c) z H yl 

x^(d 2 n<^-d^nw)=(i/cb,+d 2 )fl, 1 ~^^H yl --d x £r 2l . (25} 

The above boundary conditions give four equations in four unknowns A x , A 2 , B h and B 2 ~ 
It can be shown that this system of four equations can be reduced to just two equations in 
two unknowns as follows: 
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Once the coefficients ^ and A 2 are found, the remaining coefficients B x and i? 2 can be deter- 
mined from (25). In fact we find 
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The explicit determination of the four boundary coefficients A u A 2) B u and B 2 is now a 
straightforward although a tedious matter. The results would necessarily be lengthy and 
probably not very useful. In the next section we shall introduce a high frequency approxi- 
mation which will simplify the results a great deal without the loss of the essential aspects of 
the problem. 

2.2. High Frequency Approximation 

Any approximation that we may introduce to simplify the field components is necessarily 
contingent upon the simplification of the propagation factors in the ^-direction ; i.e., Si and s 2 
given by (12). A useful approximation will prove itself to be one based on the assumption 
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which can be brought about by either a weak magnetostatic field or by a sufficiently high wave 
frequenc} 7 . 4 Furthermore, for a high enough wave frequency it usually happens that 
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is also satisfied. If these two assumptions hold, it can be shown that the components of the 
permittivity tensor in (4) can be simplified as follows: 
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where we retained only first order terms in co^/co and v/a>. Consistent with these approxi- 
mations it can be shown that s lt2 in (12) simplifies to 

sl^s^Kk^ (32) 



4 In the earth's ionosphere, the region of validity of this approximation will be limited to frequencies in excess of 1 Mc/s. 
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where s=(kf— a 2 x —af) and ki=-^Jk . Now in regions where |s 2 |>>|k&i<*i| we can approximate 
further by writing 

nkiOLi 



Si t2 ^S±- 



2s 



(33) 



It will appear that last form of s it2 is indeed the only useful one. The usefulness of it is, however, 
reduced by the fact that it is not valid in the vicinity of s=0 which corresponds to the critical 
angle in the absence of the magnetostatic field. 

Introduction of approximation (32) into (27) and (28) results in the following: 
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(34) 
(35a) 
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(35c) 

(35d) 

(36a) 
(36b) 



As a consequence of the above results, the coefficients B x and B 2 pertaining to the air region 
can be found as follows 
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In the above expressions the first terms represent the regular field in the absence of the magneto- 
static field. The terms preceded by k represent the first order correction for the presence of 
the magnetostatic field. 

To facilitate the approximate evaluation of the various integrals we shall define certain 
fundamental integrals from which all others could be derived by differentiation. We define 
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The above integrals are not independent. Indeed, one can show that they are related as 

follows: 5 

2i 

Kq 

2id y U 2 =(d fl +^)U 3 . (39) 

In what follows we shall avail ourselves of the above relationships to check the differentiability 
of the results after integration is performed. Employing the above definitions of U u U 2 , and 
c7 3 one can write the components of the Hertzian vector in terms of them. One obtains 

+i *[i=f ( d »-^)(^+fW+ai] ^\)} (40) 

The integrals in (38) are surface integrals over the entire a { — a 2 plane of a form that 
lends itself readily to a transformation to cylindrical coordinates in the configuration space 
as well as in the transform space. Thus we employ the transformation 

x=p cos <p 

y=ps'm<p (41) 

for the space coordinates and 

a\ = X cos a 

a 2 =\ sin <r (42) 

for the transform variables. Now the fundamental integrals are of two distinct types. The 
first one is 
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The second distinct integral is of the form 
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In the appendix we show that 
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In view of the above we can rewrite the fundamental integrals as follows: 

I /»« ish 

Ui=y- -^e is o"H^(Kp)\d\ (47a) 



5 It may be shown that the differentiation under the integral sign is permissible by virtue of the fact that the resulting asymptotic expan 
sions are identical. 

509 



oo f* °° fi ish 

U 2 =k S (-1)' sin (2v+l)<p , " e<*o°Hg> +1 (\ P )d\ (47b) 

u * =k °f" a 8(8+*8 y e "° w (xp)xrfx (47c) 

where we used a well-known relationship between the Bessel and Hankel functions. 
Now w T e transform to spherical coordinates using the transformation 

z=r cos 6 

p—r sin 6 (48) 

for the configuration space and 

X=& sin (49) 

for the transformation space. The results are 

Ui= f sinPcosflff^fopsinfl) ^ vp^iv ~. (*-*> dj8 (50a) 

Jr! Vf-sin 2 )8+f cos/3 
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where ri is the appropriate path of integration in the complex /3-plane and 

Hg>(k p sin p)=H£ ) (k oP sin p) e - ik op**e. (51) 

In what follows we shall evaluate the integrals in (50) by the method of steepest descent which 
is contingent upon being able to deform the original path of integration I\ to the path of steepest 
descent I\ In the process of the continuous path deformation we must make sure that any 
singularities between these two paths are properly accounted for. To this end we shall investi- 
gate the location of the various singularities of the integrands of U u U 2 , and ?7 3 in the complex 
/3-plane. 

The examination of the integrands in (50) reveals a possibility of poles at points 
where the denominator vanishes, i.e., where sin jff=±(f)»(l + f)~*. The corresponding 
Biemann surface consists of two sheets and we shall denote the upper sheet one on which 
Im {(f — sin 2 /3)*}>0 and the lower sheet one on which Im {(f— sin 2 &)>} <0. Then it follows 
[Tyras, 1962] that when |Vf |<1 the resulting poles P 1} P 2 , P 3 , and P 4 are located in the com- 
plex /3-plane as shown in figure 2 where P 1 and P 2 are the upper sheet poles and P 3 and P 4 
are the lower sheet poles. In particular it can be shown [Tyras, 1962] that the poles P 3 and 
P 4 are on the same side of the branch cut Im{(f— sin 2 /5) * } = as the origin. 

The integrands in (50) also contain the radical (f — sin 2 ffi as a result of which the point 
j8= =fc arc sin VT is a branch point. The convergence of the integrals is assured if the path 
of integration at least begins and ends on the upper sheet. The two sheets are joined along 
the lines Im{(f— sin 2 fi)?=0 starting at the branch points. These lines are shown dashed in 
figure 2. 

We now investigate the possibility of replacing the original path of integration Ti with 
the path r corresponding to the path of the steepest descents through the saddle point 0=6. 
Since the poles P 3 and P 4 lie on the lower sheet they need not be taken into account when 
deforming the path of integration Fi into I\ The upper sheet poles P x and P 2 need not be 
taken into account either since they will never be crossed lying beyond the region of interest, 
-7r/2<0<7r/2. 
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Figure 2. The complex (3-plane. 




Next we focus our attention on the following. As we remarked earlier, the integrals 
in (50) are double valued since they contain the radical (f— sin 2 /3)*. The original path of 
integration passes over the upper sheet of the corresponding Riemann surface and it can be 
deformed into the path of the steepest descent r only when at least the beginning and the end 
lie on this sheet. In the case when the angle 6 does not exceed the critical angle 5 =arc sin Vf, 
the transition from the path Tj to the path r is accomplished without complications and the 
only contribution to the integrals will be from the saddle point. When 6 exceeds the critical 
angle 6 B a more complicated path must be devised which results in the contribution from the 
saddle point and the contribution from the integration along the borders of the branch cut. 
The latter contribution results in the phenomenon of lateral waves [Brekhovskikh, 1960; 
Tyras et al., 1963]. In the present formulation we introduced the high frequency approxima- 
tion in (33) excluding the neighborhood of the point f=sin 2 f$ from the region of validity. 
This makes the integral expressions of (50) not suitable for carrying out the branch cut inte- 
gration and thus the evaluation of the lateral wave field. 6 

The evaluation of the integrals in (50) by the method of saddle-point integration is a 
straightforward process [Brekhovskikh, I960]. We obtain for the leading term 



TJv 



J2 COS fl^MVf-sin 2 9 e ik * r 

Vf— sin 2 + f cos 6 k r 



(52) 



U, 



?"2 cos 0e lk& 



ifc /iyV-sin 2 & 



At 



Vf — sin 2 0(Vf— sin 2 0+f cos 0) 2 k r 



(53) 



The evaluation of the integral U 2 in (50b) is more complicated and we shall use the saddle- 
point method for double integration to accomplish it [Bafios and Wesley, 1953]. First we put 



F(0)e ik o* 



cos ifi-0)(ip 



where 



Now we put 



w 



cos ffff&iflfcop sin fle'Wr-rin'g 
Vf-sm 2 0(y'r-sm 2 /3+f cos 0)' 



w=p—t 



(54) 

(55) 
(56) 



6 It may be of interest to point out that in the rigorous formulation the branch points occur at si,2=0 corresponding to sin /3=±V e ± 7 ?- As 
co///o)~*0 these two branch points coalesce into one at sin (3=±ylf. The evaluation of the lateral field is possible in the case of line sources [Tyras 
et al., 1963]. 
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which transfers the saddle-point to the point w=0 in the complex w-plane. Following the 
procedure of Bafios and Wesley [1953] we make the next transformation as follows: 

2 

—=ik r (1— cos w) (57) 

and note 

<h» = x ( 58) 

dx ik r sin w 

It can be shown that w can be expanded in a power series of x as follows: 

W= (ik r)^ + 24: ' (iV) 3/2+ 640 ' (ik r)^ 2 + ' ' ' (59) 

from which we note that w is an odd function of x, i.e., w(x) = —w(—x). Using the above 
transformation we can recast the integral (54) in the form 

I=e^r G{x)e-* l2 dx (60) 

where 

v v dw 

G(x) = [F(w+d)Hl 1) (k oP sin (w+d^+Fi-w+d)!!^^ sin (-w+d))} ~ (61) 

Putting this expression under the integral sign in (60) results in two integrals with the inte- 
gration range from to +°°. In the second integral we make the substitution x=— x' and 
reverse the limits of integration. Due to the fact that dw/dx is an even function of x, we can 
add these two integrals to form a sinole one as follows: 



^ P 

ik r J - 



xF(w+e) #(i )( £ op sin ( w+e )) e -^'2 dXm ( 62 ) 



sin w 
Now we make use of the integral representation of the Hankel function 



gp>(g)= rJl'*'?', (" tf v (*iz-y*) v - 1,2 e-* 2 dy (63) 

Jir2PT(v+12)z v Jo 



Vtt2 3j T(v+1/2)s 
which we put in (62) to obtain 



/ 

ik { 



where 



A.pi(k r-vir) C™ r°° 9/ 9 

n Q(x,y)e-* 2/2 e-v 2/2 dxdv 64) 

• r^2 3 *I>+l/2)(M' J— Jo 

^_ xy*»F(w+e)l4ik p sin {w+e)-y 2 ]^ 2 . 

W x >y)- sin w [sin (w+0)]> _ * (b5) 



We recognize the integral (64) as being analogous in form to single integral of the type 
to which W^atson's Lemma can be applied. A theoretical basis extending Watson's Lemma to 
double integrals has been provided by Bafios and Wesley [1953]. Thus, at this point all that 
remains to be done to find the saddle-point contribution is to expand the integrand Q(x,y) in 
a double power series and integrate term by term. When performing such an integration we 
avail ourselves of the following w T ell-known result [Brekhovskikh, 1960] 

p &'e-**'*dx=T (*+D (IT 1 " 172 (66) 

Combining the above results we obtain for the leading term 



TT -2 cos 0e lk <M' sin V** r JU . , , -. , A7 x 

U 2 ^ , . , r • v- — V7 - ' 2 sin (2v + l)<p. (67) 

sin ^Vf — sm 2 0(vf — sin 2 0+f cos d)k r o 
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The series on the right in (67) can be summed. We write formally 

S sin (2v+l)cp=Tm (««'S« flpf Vin- 
where we used formula for summation of a binomial series. Thus 

— cos0e ifc °*^ f ~ sin2 ' e ik ° T 



U 2 



sin sin <p^— sin 2 (Vf — sin 2 0+f cos (?) ^ 



(68) 



(69) 



Collecting the above results we find the components of the Hertzian vector and hence the 
field components. Thus 



Ee, 



ikesm cos Be 



ik h^{- 



27r(f cos0+Vf— sin 2 0) 



{[" sin 2 <p 

t sin ^+^ sin ~ — , 
Lf cos 0+vr-si 



sin 2 



cos 2 <p COS 



i&o^ cos : 



2-Vf— sin 2 0(cos 0+Vf-sin 2 0) 



vyyz 



(70) 



#, 



'{' 



-ik m cos v? cos Q e ik o h M^sm 2 

V— t — < Vf — sm 2 0(cos0 Vf — sm 2 0+sm 2 6) 

27r(fcos0+Vf-sin 2 0) 



. . sin sin <p Tsin 2 0+cos dJt: — sin 2 



:[ ! 



Vf — sin 2 6> L f cos 0+Vf — sin 2 2(cos 6»+Vf-sin 2 e) 

^(f cos 0+ Vf — sin 2 e) "|\ e** r 
2(cos0+Vf-sin 2 0) J/ r 



# r n~0. 



(71) 



(72) 



One observes that the above field components contain the regular field that would exist 
in absence of the magnetostatic field and a first order correction for the presence of the mag- 
netostatic field (the terms preceded by tj). It will be seen that these corrections are not in 
general symmetric about the 0=0 and <p=0 axes. The time average Poynting vector given by 



Sr0 =2Z (l^o| 2 +l#,o| 2 )] 



(73) 



where Z is the free space impedance, is plotted in figures 3 and 4 when <p=0 and x/2 respec- 
tively for a plasma roughly equivalent to that in the lower edge of the ionosphere and the 
magnetostatic field equivalent to that of the earth. It will be noted that the correction due 
to the magnetostatic field is noticeable only in the plane <p=0 and within the region of validity 
of (29) it is significant only in the proximity of the interface. 



Figure 3. Power pattern in air of a horizontal 
magnetic dipole in magnetoplasma; <p=0°, N = 750 
electrons per cubic centimeter, h = 1000 meters, 
i> = 10 6 collisions per second. 
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Figure 4. Power pattern in air of a horizontal 
magnetic dipole in magnetoplasma; <p = irl2 f N = 7<50 
electrons per cubic centimeter, \\=1000 meters, 
v — 10 8 collisions per second. 
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3. Dipole in Air 
3.1. Rigorous Formulation 

Since this problem is in many respects similar to the preceding one, we shall not go into 
it in great detail. 

The geometry of the present problem is shown in figure 5. As before we shall work with 
the Hertzian vector of the magnetic type in the air region and with the actual field compo- 
nents in the magnetoplasma region. 



For the air region we can write immediately 



m M > 



(74) 



oo 

-'« 8irVo J J 



(75) 



whereas for the magnetoplasma region the integral representation of the field components is 
identical to that in (23). Application of the boundary conditions (25) gives a system of equa- 
tions as follows: 

U G22 ^23 ^24 

^32 C/33 C34 

1 634 6 44> 



where 



B 2 





(i 



C Vi 



_ a2[$i(gi)+ai$2Csi)] + ^X^iai 
a2$i(Si) 



(76) 



(77a) 



CH=Cl3(Sl=S2) 



(77b) 



622 



*0 



(77c) 



c = ct 2 [$i(s l ) — {kl-a\)<$> 2 (s l )} — JKX{kl(k 2 — o$Si 
a 2 $>i(Si) 



624 — ^23(^1 — #2) 



(77d) 
(77e) 
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FIELD POINT 




Figure 5. Geometry of the problem of magnetic 
dipole above a magnetoplasma half space. 



FIELD POINT 




Figure 6. Geometry of the source and the image. 



< 32 — — 
CLi 



(77f) 



C« 



a 2 #i($i) 



(77g) 



C Zi =C 33 (s 1 =s 2 ) 



(77h) 



r = ^a2^i(si)(s 1 — iKa 2 )+a\(Xs l -iKa 2 ^[a2^2(s l )+iKX^k 2 Q s l ] a\(x— f)[.<?i<i>i(si) — i^kla 2 ] nr . 



Cu=Cu(s 1 =s 2 ) 



(77j) 



b=- 



So 



(77k) 



The coefficients A lt2 and B l2 are thus determinable from (76) but as before we shall forego 
finding their explicit form in the present formulation. 

3.2. High Frequency Approximation 
Using the approximations (33) we can find the coefficients B it2 as follows: 



r 1 2f r 2a 2 

5l= \-^ + m^> +/r7 L(HFf% 



=-{ 



<4s a 



B 



2(f-l) ai 

(s+s )(s+fs ) 



+ i 



So)(«+fs )J/ 

tti(gg+a|)1 ~]\ 
V L*(«+fSo) 2 a48(* + «o)(*+f«o)J / ' 



(78) 



(79) 



To evaluate the field components we avail ourselves of the integrals U u U 2 , and U z in (38) and 
the results in (52), (53), and (69) and obtain 

515 



E$t 



-imk f sin <p /e ik ° r e ik ^ r ° \ 
2tt 1 ~J~ \~r n~J 



f cos sin <p e i1c * T 



\ 2 \ r r J ' f cos 0+Vf — sin 2 r 

^ sin cos sin 2 <p cos cos 2 <p 



h 



E, 



pO 



f cos 0+Vf— sin 2 Lf cos 0+Vf — sin 2 2Vf — sin 2 (cos 0+ V? — sin 2 0) 

ik m cos <p cos / 1 /V*o r g ft o r o \ Vf — sin 2 (sin 2 0+cos V f — sin 2 0) g**o 
2x \2V"r ~V/ + ~ 



]?} 



(80) 



f cos 0+Vf — sin 2 



irj sin 9 sin <p 



fsin 2 0+cos ^Vf — sin 2 



Vf— sin 2 e(f cos 0+Vf — sin 2 0) |_ f cos 0+Vf— sin 2 2(cos 0+Vf— sin 2 0) 



•]£} 



s r0 ~o 



(81) 



(82) 



where r and r are the distances from the source and its image to the field point respectively 
as shown in figure 6. 



4. Appendix 



The integral under consideration is 



'<, 



gi\p cos (a— <p) 

sin <t 



da. 



(Al) 



Expanding the numerator in terms of Bessel functions [Stratton, 1941] and interchanging the 
order of summation and integration gives 



<p f*ic pin<T 

J=S i n e- in *J n &p) -? da. 

±Z J-.S1IKJ 



(A2) 



Now we focus our attention on the integral on the right. By making the transformation 
z=e u one obtains 



\=\ 4 dc=2\ 

J_ T sin(r J c 



/i= 



z n dz 



sin <r 
where C is a unit circle. Now we consider 



(z+l)(2-l) 



K 



-i 



z n dz 



)r{z+l)(z-l) 
where the contour is shown in figure 7. Integrating formally it can be shown that 



(A3) 



(A4) 




Figure 7. The complex z-pla?ie. 
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/!=-2 lim | f +f } (A5) 

«->o L J r 2 J r 4 J 

which gives 

7 1 =t*[1-(-1)«]. (A6) 

Substituting the (A6) into (A2) we obtain the desired result 



/: 



t s>i\p cos (<r—<f>) 



d<r=4wi S (- 1 )" si n (2n+l)<pJ 2n+l (\p). (A7) 

sin a - o 
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